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Magnetic topology and current channels in plasmas with 
toroidal current density inversions 

D. Circfl and I. L. CaldasQ 

Departamento de Fisica Aplicada, Universidade de Sao Paulo, 05508-090, Sao Paulo, Brazil. 

The equilibrium magnetic field inside axisymmetric plasmas with inversions on the toroidal cur- 
rent density is studied. Structurally stable non-nested magnetic surfaces are considered. For any 
inversion in the internal current density the magnetic families define several positive current chan- 
nels about a central negative one. A general expression relating the positive and negative currents 
is derived in terms of a topological anisotropy parameter. Next, an analytical local solution for 
the poloidal magnetic flux is derived and shown compatible with current hollow magnetic pitch 
measurements shown in the literature. Finally, the analytical solution exhibits non-nested magnetic 
families with positive anisotropy, indicating that the current inside the positive channels have at 
least twice the magnitude of the central one. 

PACS numbers: 52.55.-s, 52.30.Cv, 41.20.Gz 



In tokamak devices a toroidal magnetic field confine 
the orbits of the charged particles inside a chamber. A 
toroidal current flowing within the plasma produces a 
complementary magnetic field preventing particle drift 
looses. Non-inductive drive mechanisms help to sustain 
this current for long pulses with small or negative induc- 
tive drive and during a slow transition from positive to 
negative toroidal current [1,] . A relevant question arising 
from this situations is that of the structure of the mag- 
netic field if the toroidal current density becomes neg- 
ative in some region of the plasma. In the last decade, 
the achievement of quasi-steady-state alternating current 
scenarios [2| and the observation of stiff structures with 
nearly zero magnetic pitch angles in a finite region about 
the plasma center [3, |j] has attracted attention to the 
problem of current reversal equilibrium configurations 
(CRECs). 

In this work we study the equilibrium topology of the 
magnetic field subjected to azimuthal (toroidal) current 
density inversions in an axisymmetric plasma. In the fol- 
lowing pages it is shown that the equilibrium topology is 
composed of non-nested families of nested magnetic sur- 
faces, where each magnetic family defines a current chan- 
nel inside the plasma. The relation between the currents 
in the channels is studied in terms of topological quan- 
tities. In a more quantitative approach, an analytical 
solution of the equilibrium problem about a region of in- 
terest provides the topology of the magnetic surfaces and 
the relevant control parameters as well as their bifurca- 
tion values controlling the transition between different 
equilibrium configurations. This is done without speci- 
fying further plasma profiles or arbitrary functions. The 
obtained solution agrees with several published equilib- 
ria [5|-l8| while providing a simple understanding of the 
control parameters and their relation with the current in 
the different channels. 

For axisymmetric systems the equilibrium magnetic 



field may be written 



B = W x V^ + FV(j>, 



(1) 



where i/j(R,z) = RA^R, z) and F(R,z) = RB^{R,z) 
are proportional to the poloidal magnetic flux and cur- 
rent respectively. B$(R,z) and A$(R,z) are the az- 
imuthal components of the magnetic field and vector 
potential and (i?, 4>, z) the usual cylindrical coordinates. 
The term "Vip x V(f> corresponds to the poloidal magnetic 
field B pi i.e. the projection of the magnetic field over a 
plane <j> = const. It is clear that B ■ V-0 = 0, so that the 
magnetic field lines remain attached to the level sets of 
ip(R, z) called magnetic surfaces. 

For axisymmetric magnetic confinement devices the 
main source of poloidal field is the plasma current, con- 
sequently, the level sets of i/){R, z) form families of nested 
tori inside the plasma. In equilibrium, the toroidal cur- 
rent flowing inside a magnetic surface labeled by ip is 
given by the Ampere's law as 

fM)It(ip)=f B-dl = (j>-f Vij)xdL = ±j> \Vip\dl, 

Jr^p j^tp -^r^, 

(2) 
where T^ is a magnetic circuit resulting from the inter- 
section of the magnetic surface ip with an arbitrary plane 
<f> = const. For a simple magnetic circuit the product 
V"0 x dl always points in the same direction, the current 
is negative when V?/> is inwards and positive otherwise. 
The existence of current density reversals suggest that 
some magnetic surface ip® contains a vanishing current; 
from (J5J this requires |V-0| = along the continuous 
magnetic circuit Tq = T^, . This requires that for any co- 
ordinate system {u, v}, the equation d u ip(u, v) = leads 
to the same relation between u and v that the equation 
d v ip(u,v) = 0. In addition, such relation must satisfy 
ip(u,v) — tpa- This kind of degeneracy is possible for 
one-dimensional problems, but leads to structural insta- 
bility in general two-dimensional equilibria [7|, i.e. any 
variation in the poloidal field, will destroy the topology 



of the surface ip = tpo. 

A more feasible situation, i.e. structurally stable, re- 
quires that VV> vanishes at isolated points inside the 
plasma. Some of these points must be saddle, each one 
introducing four branches where the condition ift(R, z) — 
ip c is satisfied, with tp c the value of the poloidal flux on 
the critical point. This guarantees that the poloidal field 
reverses over the surface ip(R, z) = ip c after vanishing at 
each saddle. 

The CRECs topologies may be classified depending on 
how the separatrix r c = Tv, , connects two branches of 
the same hyperbolic point [9] . If T c connects two oppo- 
site branches of a saddle the resulting circuit is simple 
but non-unique, i.e. there is a second simple circuit T' c 
connecting the remaining two branches; this leads to an 
even number of axisymmetric islands (Fig. QJ,). On the 
other hand, when T c connects two non-opposite branches 
the circuit is non-simple but unique and leads to an odd 
number of islands (Fig. [TJj). Separatrixes delimit sev- 
eral families of nested magnetic surfaces acting as current 
channels inside the plasma. These axisymmetric chan- 
nels arc also called non-nested surfaces in the literature. 
In Fig. Q] the poloidal field direction reveals the sign of 
the toroidal current inside each channel. In summary, 
any internal toroidal current density reversal in a two- 
dimensional magnetic equilibrium leads to a chain of ax- 
isymmetric islands with positive current about a central 
structure with negative current. 



simple circuits -{Ti,^} (Fig. [2]d) it can be shown that 





FIG. 1. Odd and even systems of axisymmetric islands. Cor- 
responding diamonds are connected by the separatrix F c . The 
arrows indicate the direction of B p and " ■" or " x" in the mag- 
netic axes indicate the current direction for each channel. 



As the poloidal field only vanishes at isolated points 
and never reverses in a regular magnetic surface, its line 
integral on a magnetic circuit can never be zero. From 
this, the current enclosed by a magnetic circuit just out- 
side the positive channels is positive, in consequence, the 
total current in the positive channels must exceed the 
current of the negative one, not only equate its magni- 
tude. Each positive channel contains an elliptic point and 
is bounded by a separatrix and two hyperbolic points [lOj ■ 
We can always build a curve 7 parallel to V-0 that passes 
through all the hyperbolic and elliptic points defining the 
positive current channels (Fig. [2^,). The circuit 7 encloses 
a vanishing current, given that B p -dl = in all its points. 
By decomposing the critic magnetic circuit r c into two 



/ B p di = J2 r x = h , i B p di = Y^r 2 = i 2 , 



(3) 



where I\ 2 is the current flowing through the region be- 
tween 7 and T12 on the i'th positive current channel. 
The relative difference between the half currents on the 
i'th channel is 7?; = (I| - I\)/I\ with /* = /» + I| the 
total current inside that channel. The current flowing 
inside all the positive channels is I + = I\ + I2 and the 
one flowing in the central channel is /_ = —I\. Using 
the previous relations we have 



F 



1-Tj' 



(4) 



with 77 = YliVi^i/Y^i^i- Equation (fj| relates the total 
current in the positive channels with the negative cur- 
rent in the central one through the parameter 77, that 
measures the anisotropy of the positive current channels. 
For small anisotropies we have i+ ~ 2J_ and in simple 
cases where the islands are created in regions with mono- 
tonic variation of the current density we expect r\ > 0, 
leading to 1+ > 2|J_|. 



(a) 





FIG. 2. (a) Regions for the half- currents between 7 and the 
separatrix. (b) Decomposition of a self-intersecting circuit 
into the simple circuits Ti and F2. Open arrows show the cir- 
cuits orientations and simple ones the poloidal field direction. 

To study the properties of the equilibrium solutions, 
like the size or number of islands, we put n in terms 
of geometric parameters. For this, we develop a local 
successive approximations method that account for the 
non-nested topologies while keeping a simplified physical 
picture of the equilibrium. 

From ((T|) and Vx£ = fi j it can be verified that 



RV ■ {R- 2 Vip) 



-MOJ0- 



(5) 



Also, from the single-fluid MHD equilibrium [ll| the force 
balance condition Vp = j x B, leads to the relation 



IMiU = MoiVWO + R^Fi^F'Oj), 



(6) 



where the prime denotes d/dip, p(ip) is the kinetic 
pressure and F(ip) is defined after (pj. Equating (0 
and ((6]) give s the well known Grad-Shafranov (G-S) 
equation 12|, |l3(, a nonlinear elliptic partial differen- 
tial equation. In ©, self-consistent surface functions 
{p(ip) , F {ip)} and boundary conditions characterize the 



magnetohydrodynamic equilibrium reached by the sys- 
tem 11] . In regular situations the equilibrium is de- 
scribed by a single family of magnetic surfaces, requir- 
ing a single choice of {p(?p),F (i/OI an d a single bound- 
ary (the plasma edge). This may not be appropriate for 
CRECs, but in a first approach a single choice of the 
sources p(tp) and F(ij)) can present current density inver- 
sions and multiple magnetic families [5|, l6j . Other works 
uses successive approximations to the solution with pre- 
scribed zero-order current density models and boundary 
conditions @,ll]. 

Formally, a broad range of choices for the arbitrary 
functions p(tp) and F(ip) or zero-order profiles may lead 
to current reversals and non-nested configurations, but 
the underlying description of the equilibrium topology 
is not restricted to the particularly chosen model. In 
the following we assume directly the existence of a small 
negative minimum of the current density in some small 
region of the plasma. With this approach we cover a wide 
range of arbitrary choices without specifying the global 
configuration of the plasma. 

Near a current density minimum at ro, a usual Taylor 
expansion gives j^r) « j^ro) + l/2(8r ■ V)l„j<p(r) + 
0(Sr 3 ). For an up-down symmetric equilibrium this can 
be casted like j^ = jo + u{r 2 — nr 2 cos 2 9), where (r, 6) 
are local polar coordinates centered at ro. Here, jo — 
j(r ) < 0, r = |i — r | and 9 is measured counterclockwise 
between en and r — Tq. The parameters i and k are 
related to the curvature and ellipticity of j(r) about r . 
For k — 0, a circle with radius a = \f— 2j /l contains a 
zero current and for r < a/y/2 we have j^ < 0. Since we 
are interested in the topology of the magnetic surfaces in 
a region containing an almost vanishing toroidal current 
we define a region of interest with length 2a, centered at 
ro (Fig®. 




Q region of interest 
( ) plasma 

I U < o 




FIG. 3. Representation of the region of interest. 

Provided that the region of interest is small we can 
solve ([5]) by a local scheme analogous to the successive 
approximations method 11]. In our case the inverse as- 



pect ratio, e = a/Ro, is a reliably small parameter defined 
from the region of interest instead of the plasma radius, 
in consequence, a fast convergence of the approximations 
is expected. Since the local current description is not in- 
tended to provide the vanishing of the current density at 
the plasma edge, the following solution is valid for r < a. 



Defining R = R + ax and z 
about the current minimum like 



<> 



d y-T 



ex 



ay we can write (O 



9 I U = (l+et)(l-2r i +2/t4 (7) 



with r measured in units of a and ip in units of fj,o\jo\a 2 Ro- 
The parameter k corresponds to the squared eccentricity 
of the current density level sets. The solutions of ([7} 
will depend on two parameters only, allowing a detailed 
study of the bifurcations that change the topology of the 
equilibrium in the region of interest. 

Assuming a large aspect ratio e << 1 and small el- 
lipticity |k| << 1, we write the nondimensional flux as 
il){r,0) = 4>o(r) + eipi(r,x) + 0(e 2 ). Excluding the el- 
lipticity in the zero-order calculations we have ipo( r ) — 
(1 — r 2 /2)r 2 /4, consequently, the first order problem be- 
comes 



AV>i 



3 - hr 2 



2 k 



(8) 



with A the 
ing Tpi(r,x) 



Laplace operator for a plane. Dcfin- 
= tpi(r,x) + nx^/Ge we have Ai/ji = 
(3 — 5r 2 ) x/2. Then, we reduce the problem to a sim- 
ple ODE by introducing the ansatz ipi = fi r )x- After 
straightforward integrations and substitutions we obtain 
the poloidal flux function to a first order in e. 



rb{r,x) = 



1 - -r 2 ) r — + ex [ 1 - % 2 



16" + 6* ■ (9) 



With this local solution we can study the magnetic field 
in the region of interest r < 0(1). For instance, if wc 
write the toroidal field as B t — Bq/(1 + ex), the mag- 
netic pitch angle, defined by 7 m = arctan(_B p /_B t ), be- 
comes 7 m = arctan(r|V?/'|), with r = /io a |Jo|/-E>o- By 
a suitable choice of the region of interest and adjust of 
k we can fit experimental measurements of j m about a 
current minimum and model the corresponding magnetic 
surfaces (Fig. [4]). This local solution agrees well with the 
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FIG. 4. Reproduction of experimental magnetic pitch angles 
in a current hole and reconstruction of the magnetic surfaces. 
Experimental data is from N.C. Hawkees et at, Phys. Rev. 
Lett. 87, 115001 (2001). 

reported Motional Stark Effect measurements in a cur- 
rent hole in [3j and predicts jo ~ —42.6 kA/m 2 , a value 
within the error reported by Hawkees et al. 



The structure of the magnetic surfaces depend on the 
relative positions of the critical points (VV> = 0), their 
type (elliptic or hyperbolic) and the value of ip on them. 
Since r and x are independent coordinates for ^ {0, 7r}, 
the off-axis critical points satisfy d r ip = d x ip = 0. The 
origin r — 0, is always an elliptic point between two crit- 
ical points in the x-axis with x c ~ ±1. The type of 
the these points depend on the ellipticity K. To under- 
stand this, we use the equation resulting of d r ip — and 
the condition for off-axis critical points \x c \ < r c . Tak- 
ing into account the domain of the first term in (0) we 
obtain |k| > e/8. The increase of the ellipticity leads 
to bifurcations about |«| ~ e/8, where off-axis critical 
points are created from those in the x-axis that change 
of type after the bifurcation occurs. A negative (positive) 
ellipticity represents a vertically (horizontally) elongated 
current density. For k < e/8, there is a single positive 
current channel about the negative one (Fig. [SJi-c). For 
K < ~ e /8 the positive channel contains an internal sep- 
aratrix (Fig. [5k,b) that for large anisotropies (or small 
e) k << —e/8 tend to split the channel by merging with 
the main separatrix (Fig. (5^,). For —e/8 < n < e/8 no 
off-axis critical points exists (Fig. |5fc) and for k > e/8, 
the on- axis hyperbolic point changes to an elliptic one 
creating a second positive channel (Fig. [5ji). 
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r\ does not correspond to a vanishing ellipticity. This is 
due to the toroidicity (represented by e), that adds an 
implicit anisotropy to the system. For this local solu- 
tion the anisotropy only vanishes for the cylindrical case 
(e = 0, k = 0) and never becomes negative, so it verifies 
i+ > 21- . The value of r\ saturates during the bifurca- 
tion at k — e/8 where, briefly, a zero current island is 
created. If we consider a very small region of interest 
(e — > 0), the current ellipticity dominates and the equi- 
librium always presents two islands. This means that, 
if a quasi-stationary transition to an internal reversed 
current density is possible, the magnetic axis is simulta- 
neously split into three separated axes. 




FIG. 5. Levels sets of i[> from Q with e = 0.1. The continuous 
separatrix defines the current channels and the dashed one the 
island internal structure. 



FIG. 6. Change in the anisotropy r\ as a function of the el- 
lipticity k for different values of e (left). Intrinsic anisotropy 
due to the toroidicity with zero ellipticity (right). 

In summary, the magnetic topology related to current 
density reversals defines several current channels within 
the plasma. The ratio between the current in the positive 
channels and the central negative current depends on a 
topological parameter measuring the anisotropy of the 
positive channels. In general terms the positive current is 
about twice the size of the central negative current. This 
anisotropy was shown to be related to the geometrical 
properties of the equilibrium in a region of interest inside 
the plasma. Then, it was observed that the anisotropy 
is always positive, indicating that the positive channels 
have at least twice the current magnitude of the negative 
one, causing the screening of this channel and forming a 
structure with net positive current. Finally, it was also 
shown that experimental magnetic pitch measures in a 
current hole are compatible with the existence of a small 
negative value of the current density. 

This work was partially supported by Conselho Na- 
tional de Desenvolvimento Cientifico e Tecnologico and 
Fundagao de Amparo a Pesquisa do Estado de Sao Paulo. 



Since the ellipticity k, and inverse aspect ratio e, de- 
fine the island sizes and the equilibrium topology, they 
must be related to the anisotropy r], in a continuous way. 
This is depicted in Fig. [6] Increasing the vertical elon- 
gation of the current density leads to a growth of the 
anisotropy of the current channels. The same is not 
true for the horizontal elongation, since the minimum 



* davidcirotaborda@gmail.com 

t ibere@if.usp.br 

[1] B. J. D. Tubbing, N. A. C. Gottardi, B. J. Green, J. A. 

How, M. Huart, R. Konig, C. G. Lowry, P. J. Lomas, 

P. Noll, J. J. O'Rourke, P. H. Rebut, D. Stork, A. Tanga, 

A. Taroni, and D. J. Ward, Nucl. Fusion 32, 967 (1992) 



[2] J. Li, J. Luo, S. Wang, P. Fu, B. Shen, F. Liu, B. Wan, [6] 

J. Shan, G. Xu, J. H. J. Yu, J. Hu, Q. Yuan, and Y. Hu, [7] 

Nucl. Fusion 47, 1071 (2007) 

[3] N. C. Hawkes, B. C. Stratton, T. Tala, C. D. Chal- [8] 

lis, G. Conway, R. DeAngelis, C. Giroud, J. Hobirk, 
E. Joffrin, P. Lomas, P. Lotte, J. Mailloux, D. Mazon, [9] 

E. Rachlew, S. Reyes-Cortes, E. Solano, and K.-D. Zas- 
trow, Phys. Rev. Lett. 87, 115001 (2001) [10] 

[4] T. Fujita, T. Suzuki, T. Oikawa, A. Isayama, T. Hatae, 

O. Naito, Y. Sakamoto, N. Hayashi, K. Hamamatsu, [11] 
S. Ide, , and H. Takenaga, Phys. Rev. Lett. 95, 075001 [12] 

(2005) [13] 

[5] A. A. Martynov, S. Medvedev, and L. Villard, Phys. Rev. 
Lett. 91, 085004 (2003) 



S. Wang, Phys. Rev. Lett. 93, 155007 (2004) 
P. Rodriguez and J. Bizarro, Phys. Rev. Lett. 95, 015001 
(2005) 

C. G. L. Martins, M. Roberto, I. L. Caldas, and F. L. 
Braga, Phys. Plasmas 18, 082508 (2011) 
The separatrix must be followed so that the trajectory is 
differentiable 

This may be generalized to include several elliptic points 
and some internal separatrix 
J. P. Freidberg, Rev. Mod. Phys. 54, 801 (1982) 
V. D. Shafranov, Sov. Phys. -JETP 26, 682 (1960) 
H. Grad and H. Rubin, in Proceedings of the Second 
United Nations International Conference on the Peace- 
ful Uses of Atomic Energy, Vol. 31 (Geneva, 1958) p. 190 



